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FORTRAN PROGRAMS FOR THE DESIGN OF LIQUID-TO-LIQUID JET PUMPS 

by Nelson L. Sanger 
Lewis Research Center 


SUMMARY 

The one -dimensional equations describing noncavitating and cavitating flow in 
liquid-to-liquid jet pumps were programmed for computer use. Each of five programs 
were written to incorporate a different set of design input conditions. The programs 
may be used for any liquid for which the physical properties are known. Calculations 
for noncavitating and cavitating performance were combined, permitting calculation of 
cavitation limits within the program. Design charts may therefore easily be developed 
without the manual iteration which is common to existing design methods. 

The equations and method of calculation are presented for each program. And in 
each case, a sample design problem is solved which illustrates the procedures and the 
types of charts that can be developed. 

The program inputs consist of pertinent pressure, flow, and geometric variables; 
estimated friction loss coefficients; and fluid properties. Outputs consist of the basic 
jet pump nondimensional parameters; other pertinent pressure, flow, and geometric 
variables; and an indication of whether the flow is cavitating or noncavitating. 

Listings of the FORTRAN IV programs are included. Execution times for each pro- 
gram are less than 1 minute on IBM-7094 equipment. 


INTRODUCTION 

The liquid-to-liquid jet pump has found increasingly wide application in recent years. 
Some examples of its diverse usage include reactor coolant circulation pumps, aircraft 
fuel pumps, and condensate boost pumps for Rankine cycle space electric power systems. 

To keep pace with the renewed interest in jet pumps, analytical and ekperimental re- 
search of their performance characteristics has also expanded. Attention has been di- 
rected toward optimization of geometry (refs. 1 and 2), cavitation performance (ref. 3), 
staged operation (ref. 4), and the operating characteristics of low-area -ratio jet pumps 


(ref. 5). Analytical and empirical relations have been developed which accurately predict 
both noncavitating and cavitating jet pump performance (refs. 3, and 6 to 10). 

Yet, despite the greater amount of information, the designer of a jet pump for a spe- 
cific application is still faced with a cumbersome task. Design charts of a general nature 
are available in some papers, but are restricted to noncavitating operation and to a rela- 
tively narrow range of area ratios. Separate calculations are necessary to check for cav- 
itation limits. And, in most cases, several manual iterations are necessary. 

To simplify and reduce the amount of work involved in the design procedure, the non- 
cavitating and cavitating procedures have been combined and programmed for computer 
use. Five design routines are presented in this report. Each of them corresponds to a 
commonly encountered jet pump design problem. FORTRAN IV listings for each are in- 
cluded. The program can be used for any liquid for which the physical properties are 
known. Therefore, for a given set of input conditions, a designer can easily and quickly 
develop a complete set of predicted performance curves showing the cavitation limits as 
well as the required physical dimensions. 


DESIGN EQUATIONS 

A schematic representation of a jet pump is shown in figure 1 , and all symbols used 
are defined in appendix A. The primary fluid (fig, 1) is pressurized by an independent 
source and leaves the nozzle as a core of high-velocity fluid. It is separated from the 
secondary stream by a region of high shear. Turbulent mixing between the two fluids 
occurs in this region, which grows in thickness with increasing axial distance from the 
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nozzle exit. The lowest pressures in the flow field occur in the shear region, and there- 
fore cavitation inception occurs there also. 


Assumptions 

The assumptions that are used in the analysis are 

(1) Both the primary and secondary fluids are incompressible. 

(2) The temperatures of the primary and secondary fluids are equal; therefore the 
specific weights are equal. 

(3) Spacing of the nozzle exit from the throat entrance is zero. 

(4) Nozzle wall thickness is zero. 

(5) Mixing is complete at the throat exit. 


Basic Parameters and Design Equations 

Four basic jet pump parameters, all expressed in dimensionless form, are used. 
They are 

(1) Nozzle-to-throat area ratio 


R 



(2) Secondary -to-primary flow ratio 


M 



(3) Head ratio 


N = 



(4) Efficiency 


r? = MN 


( 1 ) 


( 2 ) 


(3) 


(4) 
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The noncavitation analysis consists of an application of continuity, momentum, and energy 
equations across the jet pump (see ref. 8 for complete development). Because the anal- 
ysis is one -dimensional and the mixing process is three-dimensional, the analysis must 
be supplemented by empirical information to determine optimum throat lengths, nozzle 
positions, diffuser geometry, and area ratios for specific applications (e.g. , see "Design 
Considerations" section of ref. 5). 

The formula for head ratio which results from the analysis is 

2R + 2r2] ^ 2 - (1 + K. + K.)R 2 (1 + M) 2 - (1 + K ) -P- M . 2 - 

1 O X a bn 

N = ^ ■ (5) 

1 + K p - 2R - y— j^- + (1 + K t + K d )R 2 (l + M) 2 


The theoretical expression for efficiency is obtained by multiplying equation (5) by M. 
The formula for primary flow rate Wj is 


w i = 


144g 


/■ 


(P t -P 2 )144 

i/M 


/< 1 + K P >' ( 1 + K S )(^) 2 


and the formula for primary nozzle exit area A n is derived from it: 


A = 
n 


144Wjg 


g c y 






(P x -P 2 )144 


y/2g 


( 6 ) 


(7) 


Friction losses are taken into account through the use of friction loss coefficients K, 
which are based on dimensionless total -pressure losses in individual components of the 
pump, such as the primary nozzle, throat, and diffuser. The friction loss coefficients 
may be determined either by estimating the values on the basis of information in the lit- 
erature (refs. 6 to 8) or by calibrating the individual components. 

Several cavitation prediction parameters have been proposed. One of them cr^ has 
been recommended for design use in a summary report on jet pump cavitation (ref. 3). 
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It was developed independently in 1968 by this author at NASA (referred to as the alternate 
cavitation parameter a in ref. 9) , and also by Hansen and Na (referred to as a in 
ref. 10). The parameter predicts conditions at the head-rise breakdown point, which is 
also the limiting flow point (not incipient cavitation) and is defined as 


a L ~ 



where V is the secondary fluid velocity at the throat entrance (fig. 1), 
s 


v 

s ys c 


w r 


( A t - A n ) 


( 8 ) 


(9) 


A value for the minimum secondary inlet pressure required to prevent cavitation can 
be calculated from equation (8), 


p y 

2REQD - HI 2 g 


144 g_ W 2 R 


y 


+ p. 


do) 


where A n and R enter the relation from equation (1). The criterion used in the com- 
puter programs to determine cavitation-limited conditions is a comparison of P 2REQD 
and the available Pg . 

The noncavitating theory predicts experimental performance quite well over a wide 
range of area ratios and flow conditions. Comparisons between theory and experimental 
performance are presented in references 6 to 8. 

Cavitation-limited flow conditions have been investigated by various researchers, 
and empirical values for a ^ have been established. These values are summarized in 
references. A conservative design value for cr L is 1.35. Well -designed secondary in 
let regions allow values of cr L from 1.0 to 1.1 to be used; cr L is an input to each pro- 
gram and may be specified by the user . In the numerical examples presented later in 
this report, cr L = 1.1 is used. 
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DESIGN PROGRAMS AND PROCEDURES 


Because of the diverse applications possible, there are several combinations of in- 
put conditions which a designer might encounter. Five combinations are presented in 
this section and the theoretical equations are developed into five design programs . The 
equations for each program are followed by a sample design problem illustrating use of 
the program. A FORTRAN IV listing of each program is given in appendix B. Execution 
time for each program is less than 1 minute on IBM-7094 equipment. 

The choice of which program to use will depend on what input information is avail- 
able, and what output is desired. Table I summarizes the input-output features of each 
program. Until the user is familiar with all the programs, table I should be used as a 
starting point. 

Program I is one of the more versatile programs. It lends itself quite well to design 
chart development since Pj and Wg are the only inputs definitely required. The 
"varying input variables" (M, R, and Pg for program I) may be selected at random to 
permit the effects of variation of each to be investigated. If specific values for each of 

TABLE I. - INPUT AND OUTPUT VARIABLES FOR EACH DESIGN PROGRAM 


[Input common to all programs: K^, K g , K^., K^, y, p y , 


Program 

I 

II 

III 

IV 

V 

I 

II 

III 

IV 

V 

Fixed input variables 

Output variables 

p i 

P 1 

P 1 


p i 

N 

N 

N 

N 

N 

W 2 

W 2 

W 2 

W 2 


V 

V 

V 

V 

V 




P 2 

P 2 


M 






P DREQD 

P D 


w i 

w i 

w i 

w i 

w i 


d t 



d t 

P 2REQD 

P 2REQD 

P 2REQD 


P 2REQD 





d n 

P D 

P D 

P D 


P D 





R 

d n 

d n 

d 

n 

d n 


Varying input variables 

d t 


d t 

d t 


M 


M 

M 

M 





W 2 

R 

R 

R 

R 





P 1 


P 2 

P 2 
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the varying independent variables are known (either initially or from the output of another 
program) , program I may be run in a straightforward manner to produce only one set of 
output. 

Program n is used when the throat diameter of a pump is known, either as a design 
constraint, or as part of an existing pump that is to be redesigned. Program HI is the 
one program that will most often be used in conjunction with some other program. It is 
used when a pump must be designed to operate quite close to the cavitation limit. Pro- 
gram III identifies the cavitation-limited pump configurations. With this information the 
designer can apply a safety margin to the appropriate parameter and recalculate the final 
design using another program (e.g. , program I). 

Program IV is used when secondary flow rate Wg and pump pressure rise P^ - P2 
are known, and when there is some flexibility in the choice of driving pressure P^ or 
flow Wj. Finally, program V is used when the jet pump geometry is completely speci- 
fied and it is desired to know the off -design performance. 

Some design problems will probably occur which were not anticipated by these five 
computer programs. However, in most cases, it should be possible to create new pro- 
grams by combining the appropriate design equations . 

The sample design problems presented after each program do not represent the full 
range of applicability of each program. For some applications, enough information will 
be available to permit a straightforward once-through design procedure. In other cases, 
it will be necessary to create design charts before arriving at a final design. 

Similarly, in the sample problems some programs are used in conjunction with 
others (programs in with I, I with H, and IV with V). How certain programs are used 
with each other, or if they are, will also depend on the specific application. The combi- 
nations used in the sample problems in this report are not suggested as the only possi- 
bilities open to a designer. As experience is gained using the programs, the potential 
relations between them will become clearer. 

Finally, a word about design compromises. In general, the ideal jet pump design 
would possess several desirable but mutually unattainable qualities. Large amounts of 
secondary flow W 2 would be pumped by a minimum of primary flow W^. This corre- 
sponds to a high flow ratio, M = Wg/Wj. The pressure supplied by an outside source 
would be kept low while jet pump pressure rise was maximized (P D - P 2 ). This corre- 
sponds to a high head ratio, N = (P D - P 2 )/(P 1 - P D ). Efficiency would be high (77 = MN) , 
and operation would be cavitation free. 

In practice, of course, compromises must be made. Achieving all these idealized 
goals concurrently would violate the laws of conservation of energy and momentum. Some 
of the compromises encountered in practice are illustrated in the sample problems. 
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Program I 


Variables . - The known and unknown variables incorporated in program I are as fol- 
lows: The known variables are 

(1) Primary fluid inlet pressure P^ 

(2) Secondary flow rate W2 

(3) Fluid properties y and p y 

(4) Friction loss coefficients K p , K g , K t , and 

(5) Cavitation parameter cr L 

(6) Secondary inlet pressure P2, flow ratio M, and area ratio R (to be selected and 
' their ranges varied) 

The variables to be calculated are 

(1) Outlet pressure P D 

(2) Required secondary inlet pressure I^REQD 

(3) Area ratio R 

(4) Nozzle diameter d n 

(5) Throat diameter d^ 

(6) Primary flow rate 

(7) Head ratio N 

(8) Efficiency 77 

Equations . - The program uses the design equations in the order presented. From 
the input information and the definition of flow ratio (eq. (2)), the primary flow rate is 
calculated 


W 2 

W- = — - (2) 

1 M 

and from it and other input information, the primary nozzle area and diameter are cal- 
culated: 


A 


n 


144W J g 

g c r 



(1 + K p ) - (1 + K g ) 




- P 2 )144 

y/2g 



( 7 ) 
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Throat area is computed from the definition of area ratio (eq. (1)) 


A 


t 


A 



and throat diameter is calculated from 


(1) 


d 


t 



Head ratio N is computed from equation (5) and is multiplied by flow ratio M to' obtain 
efficiency: 


N = f(M, R, Kp, K g , K t , K d ) 
17 = MN 


(5) 

(4) 


Outlet pressure is calculated from the definition of head ratio (eq. (3)) and input values 
for and P 2 , 


P D = 


NP 1 + P 2 
1 + N 


(3) 


Total flow rate is calculated by summing primary and secondary flow rates, 


w T = Wj + w 2 


The minimum secondary inlet pressure required for cavitation-free operation is com- 
puted and compared to the input P 2 : 


P 2REQD “ f g 


w 2 r 


144 g 
y e c A n (1 - R) 


+ Pr 


(10) 


If P 2 is greater than P 2REQD’ the flow is noncav i tatin g and a message indicating this 
is printed out. 
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Sample design problem. - A jet pump is used to circulate 1200° F liquid sodium 
through the core of a nuclear reactor. It has a throat diameter d t of 0. 930 inch and an 
area ratio R of 0.0425. An auxiliary pump supplies a weight flow Wj of 1 pound per 
second of drive fluid at a pressure of 80 psia to the jet pump, which pumps a weight 
flow Wg of 3 pounds mass per second from an inlet pressure Pg of 25 psia. A redesign 
of the reactor core requires that the jet pump produce a pressure Pp of 30 psia instead 
of the present 28. 5 psia. This will require a new jet pump design. 

For this sample problem the known variables are P^ = 80 psia; Pg = 25 psia; 

Pp = 30 psia; Wg = 3 lbm/sec; = 1 lbm/sec; y = 49.33 lbf/ft^; p y = 0.96 psia; 

K = 0.03, K g = 0. 1, = 0. 1, and Kj = 0. 1 (estimated); and = 1. 1. The variables 

to be' calculated are R, d R , and d^. 




Area ratio, R 


(b) Efficiency and required secondary inlet pressure as 
function of area ratio. 

Figure 2. - Program I sample problem. 
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Since pressure and flow requirements are completely prescribed, the design proce- 
dure is straightforward. The results are shown in figure 2 as functions of area ratio. To 
achieve an outlet pressure of 30 psia requires a pump having an area ratio of 0.08, a 
throat diameter of 0.668 inch, and a nozzle diameter of 0. 189 inch. Such a pump will 
have an efficiency of 30.2 percent. The secondary inlet pressure required to prevent 
cavitation is 5.3 psia, well under the 25 psia available. 


Program II 

Variables . - Program II uses the following information to calculate the unknown 
variables: 

(1) Primary fluid inlet pressure Pj 

(2) Secondary flow rate Wg 

(3) Throat diameter d^. 

(4) Fluid properties y and p y 

(5) Friction loss coefficients Kp, K g , K^, and 

(6) Cavitation parameter a L 

(7) Secondary inlet pressure Pg and area ratio R (to be selected and their ranges 

varied) 

The variables to be calculated are 

(1) Outlet pressure P D 

(2) Required secondary inlet pressure ^REQD 

(3) Flow ratio M 

(4) Primary flow rate 

(5) Nozzle diameter d n 

(6) Head ratio N 

(7) Efficiency 17 

Equations . - The program calculations are performed using the following equations 
and procedure in the order listed. 

Knowing the throat diameter, the throat area is calculated, and from the definition 
of area ratio (eq. (1)), the nozzle area is computed: 

A t = 0.7854 d^ 

A n = A t R (1) 
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The only flow parameter known is the secondary flow rate W 2 . If either W ^ or M 
were known, the other could be calculated from the definition of flow ratio, M = W 2 /W 1 
(eq. (2)). This equation is used in an iterative procedure in this program to determine 
both and M. A first approximation for is made by dropping the term in equa- 
tion (6) which contains M. 


W^l) = 


yA g 

[_ n & c 

144g 



- P 2 )144 

r/2g 




(first approximation for W^) 


Having this value for W^, a corresponding first approximation for flow ratio can be cal- 
culated: 


W 

M(l) = ‘L (2) 

w 1 (i) 

An iteration loop is then begun using the complete equation (6). Each succeeding W ^ is 
computed using the flow ratio M calculated in the preceding iteration. The resulting 
value of M calculated is compared to the value calculated in the preceding iteration. 
When the percentage deviation between the two is less than 0.05, the loop is completed 
and a final value for is completed: 
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W^j) = 


( 6 ) 





- P 2 )144 

r/2g 


144g 



+ K p )- (1 


K s> 


M(j - 1)R 


1 - R 




M(j) = 



( 2 ) 


AM = M(j) - M(j - 1) 

Percent deviation = x 100 
M(j) 

If the percent deviation is greater than 0.05, recalculate the loop. 
If the percent deviation is less than 0.05, set M = M(j). 


W l = 


^n^c 

144g 


/ 


(P x - P 2 )144 

y/2g 


|/(l-K p )-(l + K s )^J 


Head ratio, efficiency, and outlet pressure are then calculated: 


N = f(M, R, K p , K g , K t , K d ) 


7) = MN 

P D = (NPj + P 2 )/(l + N) 
And finally, a cavitation check is made: 


( 6 ) 


(5) 

(4) 

(3) 
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2REQD 


u T 

J_i y 

144 2g 


w 2 r 


144 g 


r g c A n d-R) 


+ P, 


( 10 ) 


If P 2 is greater than P 2REQD’ flow is none a vita ting, and a message indicating this is 
printed. 

Sample design problem . - The original jet pump in design problem I had a diameter 
d t of 0.930 inch and an area ratio R of 0.0425. Rather than build an entirely new jet 
pump, a designer may wish to remove the nozzle from the pump body and replace it with 
one having a different diameter. 

For this sample problem the known variables a n Pj = 80 psia; P 2 = 25 psia; 

P D = 30 psia; = 3 lbm/sec; y = 49.33 lbf/ft 3 ; p y = 0.96 psia; K p = 0.03, K g = 0. 1, 



(a) Outlet pressure and nozzle diameter as function of 
area ratio. 



(b) Efficiency and primary flow rate as function of area ratio. 
Figure 3. - Program II sample problem. 
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^ = 0.1, and = 0. 1 (estimated); and a L = 1. 1. The variables to be calculated are 
R, d R , and Wj. 

The design curves corresponding to this set of conditions are presented in figure 3. 
Once again the design procedure is straightforward, but a compromise is involved. To 
achieve a higher outlet pressure than the original 28. 5 psia (see problem I), using the 
same primary inlet pressure and same throat diameter, requires a greater amount of 
primary fluid and therefore a larger nozzle diameter (0.228 in. compared to 0. 192 in.). 

The results show that a jet pump having an area ratio R of 0.06 provides an outlet 
pressure of 30.0 psia. But it requires 1.41 pounds mass per second of primary flow rate 
to do it. The design decision (comparing problems I and II) is whether the reduced cost 
and simpler replacement features are worth the required extra flow rate of 0.41 pound 
mass per second and the lower efficiency (21.3 against 30.2 percent). 


Program III 

Variables . - Program El takes the following information and calculates the corre- 
sponding flow and geometric variables: 

(1) Primary fluid inlet pressure P^ 

(2) Secondary flow rate Wg 

(3) Required outlet pressure p dreqd ’ * ower for outlet pressure 

(4) Fluid properties y and p y 

(5) Friction loss coefficients Kp, K g , K^, and 

(6) Cavitation parameter cr L 

(7) Flow ratio M and area ratio R (to be selected and their ranges varied) 

The variables to be calculated are 

(1) Primary flow rate Wj 

(2) Outlet pressure P^ 

(3) Required secondary inlet pressure f^REQD 

(4) Nozzle diameter d n 

(5) Throat diameter d^. 

(6) Head ratio N 

(7) Efficiency t] 

Equations . - The calculations are performed using the following equations in the 
order listed. From the input information, values are calculated for primary flow rate 
Wj, head ratio N, and efficiency ry. 


W 


1 



( 2 ) 
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N = f(M, R, Kp, K s , K t , K d ) 


(5) 


77 = MN (4) 

An iteration loop is then begun. A first estimate for the secondary inlet pressure P 2 is 
calculated from the definition of head ratio (eq. (3)) using the lower limit for outlet pres- 
sure ( p dreQd) ^ P D‘ Nozz l e area A n is computed from equation (7) and P 2 REQD 
from equation (10). The value thus obtained for p 2 reQD is set equal to P 2 and used 
to recalculate A fi and p 2Re q e until the difference between successive iterative calcu- 
lations for p 2 reQD is within specified limits. 

The first estimate for P 2 is 

P 2 = P DREQD " N ( p i " P DREQD) 


144Wjg l/( 1+R p) ~ (1 + K s } (t^) 
g c y /(P 1 - P 2 )144 

V y/2g 


JL 


2REQD 144 2g 


144g W 2 R 


~l2 


7Z C A n (l - R) 


+ P, 


Percent deviation = 


P - P 

2REQD _2 

P 2REQD 


x 100 


(10) 


If percent deviation is greater than 0.05, set P 2 = p 2 RE q D in equation (7) for A n , and 
recalculate the values for A n and p 2RE QD - If percent deviation is less than 0.05, con- 
tinue to the next step. 

After the loop is satisfied, the outlet pressure and geometry are calculated from 
equations (3) and (1) and the formulas for d t and d n : 

p = NP 1 + P 2REQD (3) 

D 1 + N 
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( 1 ) 



Thus, beginning with a small amount of specified information, this program computes a 
jet pump configuration designed to operate at the minimum possible secondary inlet pres- 
sure, P 2 = ^2REQD* H avan g arrived at a geometric configuration and an operating point, 
the designer may then apply a safety margin to any of several variables (e.g. , P 0 , d , 
dj_, R, or W 2 ) and enter any of the other programs to compute the final pump design. 

Sample design problem . - A Rankine cycle system is to be used for generating on- 
board electric power for a spacecraft. Working fluid is liquid potassium. A condensate 
pump is to be designed, and a jet pump will be needed to act as a booster pump to in- 
crease the inlet pressure to the condensate pump to 10 psia. A pressure of up to 300 psia 
will be available from the system to drive the jet pump (P^) , and 2. 5 pounds mass per 
second of fluid (Wg) must be supplied by the pump to the system. Fluid at the inlet of the 
jet pump will be at a temperature of 1000° F. System size requirements place a con- 
straint on throat diameter, limiting it to less than 1 inch. And radiator weight limita- 
tions specify that available secondary inlet pressure P 2 will be less than 4 psia. The 
designer must determine the amount and pressure of the recirculated fluid necessary to 
drive the jet pump (Wj and Pj), the size of the jet pump (d t , d n , and R), and the inlet 
pressure of the condensate fluid required to prevent cavitation (P 2 reQd^ * 

For this sample problem the known variables are Pj = up to 300 psia; W 2 = 2. 5 lbm/ 
sec; P dre q d = 10 psia; y = 44.38 lbf/ft 3 ; p y = 1.1 psia; K p = 0.03, K g = 0.1, K t = 0.1, 
and = 0. 1 (estimate based on literature); and a L = 1.1. The variables to be calcu- 
lated are P 2 REQD’ w i> P l> R > d n > d t’ andactual p d- 

Program HI was run for values of P^ of 100, 200, and 300 psia and over a range of 
flow ratios M from 1 to 6 and area ratios R from 0.01 to 0. 10. Results are plotted in 
figure 4. 

Figures 4(a-l) to (a -3) are plots of outlet pressure and required secondary inlet 
pressure as a function of area ratio for six values of flow ratio. 

Figures 4(a-l), (a-2), and (a-3) are for primary inlet pressures of 100, 200, and 300 
psia, respectively. In figure 4(a-l) a horizontal line has been drawn, corresponding to 
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the Pj)reQD of 10 psia. It intersects each of the six flow ratio curves at a specific 
area ratio. This area ratio and the flow ratio corresponding to it are then used to locate 
points on the ^REQD -a 8 ainst -R curves so as to construct a characteristic curve for 
P D = 10 psia. 

Figure 4(b) is a plot of efficiency as a function of area ratio for each of the six flow 
ratios. Superposed on this figure are characteristic curves for P^ of 100, 200, and 
300 psia, constructed from the intersections of the curves for p dreQD ~ ^ psia of fig- 
ure 4(a). Therefore, figure 4(b) is restricted to a P D of 10 psia. 

A study of figures 4(a) and (b) will reveal the design compromises that must be made. 
Figure 4(a) shows that as Pj increases, the characteristic curves for P D = 10 psia 
shift to lower required secondary inlet pressures, a favorable trend. But figure 4(b) 
shows that increasing P^ corresponds to decreasing efficiency, an unfavorable trend. 
Figure 4(a) also shows that to achieve the desirable goal of high -flow-ratio operation 
(low Wj) means accepting a high required secondary inlet pressure, an undesirable 
trend. So, as observed earlier, it is impossible to achieve all the desirable goals con- 
currently (in this case, high flow ratio, low required secondary inlet pressure, and high 
efficiency) . 

In selecting an operating point, it should first be recognized that the required second- 
ary inlet pressure produced by the program is the minimum possible operating pressure 
and the resulting jet pump geometry is calculated based on this minimum pressure. It is 
unlikely that this set of conditions (P 2 and geometry) will ever constitute the final de- 
sign point because of the lack of operating margin with respect to cavitation. But, having 
established the minimum possible secondary inlet pressure and an acceptable throat and 
nozzle diameter, the designer can then easily determine the final design by making one 
pass through program I, for example. The input to program I is the same as for pro- 
gram in except that a value for secondary inlet pressure must also be given. When the 
minimum secondary inlet pressure from program III is known, it can be increased by an 
appropriate margin of safety and used as input for program I. 

Figures 4(a) and (b) may be reviewed to clearly illustrate the procedure. The object 
is to keep P 2 less than 4 psia and concurrently to maximize flow ratio M. Many oper- 
ating po : nts could be chosen. Three are indicated in figures 4(a) and (b), and are listed 
below. The throat diameters corresponding to these points are noted in figure 4(c), a 
plot of throat diameter against area ratio. 

(1) Pj = 100 psia; M = 3; rj = 21.7 percent; I^REQD = psia; R = 0.0466; 

W 2 = 2.5 lbm/sec; Wj = 0.83 lbm/sec (fig. 4(a-l)), (b), and (c)). 

(2) Pj = 200 psia; M = 5; rj = 17.0 percent; P 2 REQD = 3 - 75 P sia > R = °- 0215 > 

W 2 = 2.5 lbm/sec; = 0.50 lbm/sec (fig. 4(a-2), (b), and (c)). 

(3) Pj^ = 300 psia; M = 6; t] = 13.3 percent; P 2 REQD = 3.4 psia; R = 0.0139; 

W 2 = 2.5 lbm/sec; W 1 = 0.42 lbm/sec (fig. 4(a-3), (b), and (c)). 
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For all cases, Pg is less than 4 psia (fig. 4(a)) and is less than 1 inch (fig. 4(c)). 

A compromise must be made between high flow ratio and high efficiency. At M = 6 
the specified can be pumped with one-half as much primary flow as at M = 3, there- 
by reducing the size and weight of the main stage pump. Since 300 psia is available for 
use, the choice made by the author for this set of conditions is the configuration corre- 
sponding to operating point (3). In making the choice, some efficiency was conceded. 

To arrive at the final design, program I was used. No figures are presented herein, 
but the final results are given below. The secondary inlet pressure Pg was set at the 
maximum allowable 4 psia, which provides a safety margin of greater than 15 percent 
over P 2 reQD P sia - This and the other conditions cited earlier for operating 

point (3) were used as input to program I and the results were 

(1) Area ratio, R = 0.0125 

(2) Throat diameter, d t = 0. 749 in. 

(3) Nozzle diameter, d n = 0.084 in. 

(4) Primary inlet pressure, P^ = 300 psia 

(5) Secondary inlet pressure, Pg = 4 psia 

(6) Minimum operating secondary inlet pressure, ^REQD = 2 -^ P sia 

(7) Outlet pressure, P D = 10 psia 

(8) Efficiency, tj = 12. 5 percent 

(9) Primary flow rate, = 0.42 lbm/sec 


Program IV 

Variables . - The input variables for program IV are 

(1) Secondary fluid inlet pressure P 2 

(2) Outlet pressure P^ 

(3) Secondary flow rate Wg 

(4) Fluid properties y and p y 

(5) Friction loss coefficients K^, K g , K^, and 

(6) Cavitation parameter cr L 

(7) Area ratio R and flow ratio M (to be selected and their ranges varied) 
The output of the program is 

(1) Primary fluid inlet pressure Pj 

(2) Primary flow rate 

(3) Nozzle diameter d fi 

(4) Throat diameter d^. 

(5) Head ratio N 

(6) Efficiency 77 

(7) Indication of cavitating or noncavitating flow 
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Equations . - The calculations are performed in the order that the following equations 
are listed. 

The secondary flow rate Wg is known and values for flow ratio M are selected and 
varied. Primary flow rate is then calculated from equation (2) 

W 2 

W- = -2 (2) 

1 M 


Head ratio and efficiency are then calculated, 

N = f(M, R, Kp, K g , K t , K d ) 


r] = MN 

The definition of head ratio is used to compute primary inlet pressure, 


P 1 = P D + 


P D- P 2 

N 


and equation (7) is applied to calculate primary nozzle area, 


144W t g |/d + K p) - (1 + K s> (yrg) 2 

s cV /(Pj - Pg) 144 

\ y/ 2g 


(5) 

(4) 


(3) 


(7) 


Throat area is calculated from the definition for area ratio 


A t = — 
t R 


( 1 ) 


and throat and nozzle diameters are calculated from the simple area relations 




7854 
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/ 


n 


0.7854 


Finally, a cavitation check is made by computing P 2REQD and com P arin g it with the 
available secondary inlet pressure Pg, 


p L y 

2REQD - — 2g 


144 g W 2 R 


“|2 


y Sc V 1 - R > 


+ P, 


( 10 ) 


If P 2 is greater than P 2REQD’ Row is noncay it a ting and a message indicating this is 
printed out. 

Sample design problem . - A jet pump is to be designed to pump 1.25 pounds mass 
per second of JP-4 aircraft fuel at a temperature of 100° F from a pressure of 6 to 12 
psia. No more than 200 psia will be available from the system for primary fluid pres- 
sure P 1 at design-point operating conditions. It is desired to use as little primary fluid 
as is consistent with high efficiency and cavitation -free operation. The designer must 
determine the size of the jet pump and the primary inlet pressure required. 

For this sample problem the known variables are P 2 = 6 psia; P^ = 12 psia; P^ = up 
to 200 psia; W g = 1.25 lbm/sec; y = 47.6 lbf/ft 3 ; p y = 2.0 psia; K p = 0.03, K g = 0. 1, 

K t = 0.1, and ^ = 0.1 (estimated); and cr L = 1.1. The variables to be calculated are 
Pp R, d n , and d^. 

Program IV was run for values of flow ratio M ranging from 1 to 5 and for area 
ratios R from 0.01 to 0.10. The results are plotted in figures 5(a) to (c) as a function 
of area ratio. Before selecting an operating point, and thereby a geometric configura- 
tion, the specific application must be considered. There are two general classes of ap- 
plications or operating conditions. First, if no off -design operation is expected (i.e. , 
fixed point operation) , the configurations which operate very close to the cavitation limit 
and have high efficiency can be selected. For the curves shown, this would correspond, 
for example, to R = 0. 04, M = 4 and R = 0.07, M = 3 (figs. 5(a) and (b)). 

The second class of operation would correspond to the sample problem cited; that is, 
off -design operation over a range of flow is expected. In such a case the operating point 
must be chosen to allow a margin for increase in flow W 2 without causing cavitation. 

Let us take R = 0.07, M = 3 (fig. 5(a)) for an example. If such a configuration were 
to be operated over a range of flows , it would follow the vertical operating line shown in 
figure 5(a) (constant area ratio R). Following such a path, it is found that the intersec- 
tion with the M = 4 curve is in the cavitating region. The operating range of this con- 
figuration is obviously restricted. Therefore, another point should be chosen on the non- 
cavitating portion of the same flow ratio curve. It should be located at a greater distance 
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(c) Throat and nozzle diameter as function of area ratio 
and flow ratio. Outlet pressure, P D = 12 psia. 

Figure 5. - Program IV sample problem. 



from the cavitation-limit point and such that operation at a fixed area ratio (vertical line) 
would indicate cavitation-free operation at higher flow ratios M. The point M = 3, 

R = 0.025, Pj = 148 (figs. 5(a) to (c)) would be a good example. Once such a point is de- 
termined, the predicted operating range and cavitation capabilities can be checked by 
using program V (results from program IV used as input to program V). Should the op- 
erating range appear unsatisfactory, the designer can then return to figure 5(a), select 
another point and repeat the process . 

That procedure was used to determine the design point (R = 0.025) indicated in fig- 
ures 5(a) and (b). A further explanation of the procedure will be given in the sample de- 
sign problem for program V. The design point conditions and geometry corresponding 
to the point indicated in figures 5(a) to (c) are 

(1) Primary inlet pressure, P^ = 148 psia 

(2) Flow ratio, M = 3 

(3) Area ratio, R = 0.025 

(4) Efficiency, r\ = 13.2 percent 

(5) Nozzle diameter, d R = 0.099 in. 

(6) Throat diameter, d^ = 0.624 in. 


Program V 

Variables . - Program V is used when thr jet pump geometry is known and it is de- 
sired to predict off-design performance. The known variables are 

(1) Primary inlet pressure P^ 

(2) Secondary inlet pressure Pg 

(3) Area ratio R 

(4) Nozzle diameter d n 

(5) Throat diameter d^ 

(6) Fluid properties y and p v 

(7) Friction loss coefficients K^, K g , K^, and 

(8) Cavitation parameter cr^ 

(9) Flow ratio M (to be selected and the range varied) 

The variables to be calculated by the program are 

(1) Outlet pressure P D 

(2) Required secondary inlet pressure P 2REQD 

(3) Primary flow rate 

(4) Secondary flow rate Wg 

(5) Head ratio N 

(6) Efficiency rj 
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Equations. - The equations and calculation procedure are as follows: All the infor- 
mation necessary to calculate head ratio is available as input, 


N = f(M, R, K p , K g , K t , K d ) (5) 

Outlet pressure can be calculated from the previously calculated N and from input values 
for primary and secondary inlet pressures P 4 and P 2 . 



NP 1 + P 2 
1 + N 


(3) 


2 

Nozzle area is calculated from the area formula A fl = 0. 7854d n and is used in equa- 
tion (6) with input information to compute primary flow rate 


w i = 


144g 




(Pj - P 2 )144 


f/a * 


V 


r/2g 


(1 + K s ) 


fe) ! 


Secondary flow rate is calculated from the definition of flow ratio, 

w 2 - 


( 6 ) 


( 2 ) 


and the cavitation limit is checked, 


ct t 
_ L y 


2REQD - H4 2g 


1 .. W„R 

144 g 2_ 


y g c A nd - R ) 


+ P. 


(10) 


If P 2 is greater than P 2REQD’ flow is noncavitatil ^ and a message indicating this is 
printed out. 

Sample design problem . - The jet fuel pump designed in the program IV sample 
problem must be able to operate over a range of flows and pressures to meet varying 
engine requirements. The designer is therefore interested in developing a series of 
predicted performance curves for the jet pump for the conditions and flow rate range 
specified as follows: 
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The known variables are = up to 400 psia (to meet the flow range requirement); 

Pg = 6 psia; R = 0.025; d n = 0.099 in. ; = 0.624 in. ; = 0. 75 to 2. 5 lbm/sec; 

M = 3.0; r = 47.6 lbf/ft 3 ; p y = 2.0 psia; K p = 0.03, K g = 0.1, K t = 0.1, and K d = 0.1 
(estimated); and <j l - 1.1. The variables to be calculated are P, P 2REQD’ W l’ N ’ and 
7] for any point in the operating range . 

At design conditions, the primary inlet pressure required by the jet pump will be 
148 psia; the geometric configuration was sized according to this requirement. As the 
engine operating conditions vary for off -design operation, the primary inlet pressure 
available to drive the jet pump will change and the maximum available P^ will be 
400 psia. Therefore, input values to program V for P^ were arbitrarily selected at 
50, 100, 200, 300, and 400 psia. 

The performance curves are presented in figure 6. Of most direct use to the de- 
signer is figure 6(a-l) which gives developed outlet pressure Pp and the corresponding 
flow rate Wg to the system. Outlet pressure is also plotted in figure 6 (a -2) as a func- 
tion of flow ratio. In both figures an engine operating line was constructed assuming jet 
pump operation at constant flow ratio. The indicated cavitation-limiting secondary flow 
rate Wg for all operating conditions is 2.6 pounds per second (fig. 6(a-l)) because 
secondary inlet area and pressure are fixed. However, the cavitation-limiting flow ratio, 
M = Wg/Wj, varies (fig. 6(a-2)) because changes in primary inlet pressure cause changes 
in primary flow rate . 

Nondimensional jet pump parameters N and rj are plotted as a function of flow 
ratio in figure 6(b). And for convenience, the secondary inlet pressure required to pre- 
vent cavitation is plotted in figure 6(c) for a range of flow ratios and primary inlet pres- 
sures. 

In selecting a jet pump to meet program IV sample problem requirements, the in- 
formation presented in figure 6 would be used directly. The configuration finally chosen 
(the performance of which is shown in figs. 6(a) to (c)) had the following characteristics: 

R = 0.025; M = 3; = 148 psia; d n = 0.099 inch; d t = 0.624 inch; and T) = 13.2 percent. 

The specified secondary flow rate range Wg was 0.75 to 2. 5 pounds mass per sec- 
ond. Figure 6(a-l) shows that the cavitation-limiting flow rate is 2.6 pounds mass per 
second, and this configuration is therefore acceptable. 

On the other hand, a configuration which was not acceptable corresponded to the point 
on figure 5(a) where R = 0.0185, P^ = 200 psia, and M = 4.0. This was initially attrac- 
tive because of its higher flow ratio of M = 4.0. However, when data corresponding to 
this configuration were entered into program V, the results (not shown here) gave a 
of 2.3 pounds mass per second as the limiting secondary flow rate, less than the upper 
limit of the specified flow range. 
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Input and Output 

This section describes the procedures for entering input data into the program and 
the form in which output is printed. 

Input. - All input variables for each program are entered with NAMELIST declara- 
tions. The first set of data entered is the same for each program. It consists of the 
friction loss coefficients K p , K g , K t , and K d ; the specific weight of the fluid y\ the 
vapor pressure of the secondary fluid p ; and the limiting cavitation parameter o 
This group of data is identified as CARD1, and a sample is shown below 
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On IBM-7094 equipment the set of data in a NAMELIST declaration must begin with a 
$ in card column 2, followed by the group identification name (i.e. , CARD1 in columns 
3 to 8); and a $ must appear at the end of the data on the last card. 

The second group of data entered in each program is identified as CARD2 and differs 
for each program. A list of CARD2 input variables for each program follows: 

(1) Program I: PI, W2, NORS, R, NOMS, M, NOP2, P2 

(2) Program II: PI, W2, DT, NOP2, P2, NORS, R 

(3) Program III: PI, PDREQD , W2, NOMS, M, NORS, R 

(4) Program IV: NORS, R, NOMS, M, P2, PD, W2 

(5) Program V: PI, P2, DN, DT, R, NOMS, M 

The variables in a specific group may be entered in any order. An example of input for 
the program HI sample problem is given here: 
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The NAMELIST declaration permits a simple variation of parameters (e.g. , PI 
above) with a minimum of card punching and with no concern for format fields. 

Output. - Output variables are printed according to FORMAT declarations. Perti- 
nent input information is printed first. Then program -calculated information and special 
messages are printed. A sample output sheet from the program III sample problem is 
given here. 
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CONCLUDING REMARKS 


The one -dimensional equations describing noncavitating and cavitating flow in liquid- 
to-liquid jet pumps were programmed for computer use. Each of the five programs was 
written to incorporate a different set of design input conditions. The programs may be 
used for any liquid for which the physical properties are known. Calculations for non- 
cavitating and cavitating conditions were combined, permitting calculation of cavitation 
limits within the program. Design charts may therefore easily be developed without the 
manual iteration which is common to existing design procedures. 

Three types of input data are required for each program. One is composed of fluid 
properties, friction factors, and other constants. Another type is made up of certain 
fluid dynamic and geometric parameters which are specified invariant by design require- 
ments . And the third type of input is composed of parameters which may be varied to 
allow flexibility in choice of the design point. 

The programs are adaptable in use. Single-pass design-point calculations may be 
made if the design requirements are fully specified. Or, if some of the parameters are 
variable, one or more programs may be used to construct elaborate design charts. 

Program I is a versatile program which requires only two invariant input param- 
eters: primary inlet pressure P^; and secondary flow ratio Wg. Through variation of 
flow ratio M, area ratio R, and secondary inlet pressure P 2 design charts may be 
constructed which specify the other geometric and fluid dynamic parameters. Program II 
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is used when the designer knows the throat diameter of the jet pump d^., the primary in- 
let pressure P^, and the secondary flow rate Wg. 

Program III is used when the primary inlet pressure P^, the secondary flow rate 
Wg, and the minimum allowable outlet pressure P^ are specified. It calculates a con- 
figuration which is sized for operation just at the cavitation -limiting condition. This pro- 
gram may be used to establish a first-approximation design for jet pumps which must op- 
erate close to the cavitation limit. In such cases, a cavitation safety margin is applied 
to output parameters from program III, and those data are used as input to one of the 
other programs to determine the final design. 

Program IV is used when the secondary flow ratio Wg, the secondary inlet pressure 
Pg , and outlet pressure P D are known (i . e ., Wg and jet pump pressure rise) . Pro- 
gram V is used when the jet pump geometry has been specified (area ratio R, nozzle 
diameter d n , and throat diameter d^) and the primary and secondary inlet pressures P^ 
and Pg are known. The off-design performance of an existing jet pump is calculated by 
this program. 

A sample design problem was solved for each program. In some cases, it was 
shown how two programs may be used in series. FORTRAN IV listings of the programs, 
sample input data cards, and an output data listing sheet are also included. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 18, 1971, 

128-31. 
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APPENDIX A 


SYMBOLS 


FORTRAN 

Mathematical 

Definition 

variable 

symbol 


AN 

A 

n 

2 

area of primary nozzle at nozzle exit plane, in. 

AT 

A t 

2 

area of throat, in. 

DN 

d n 

diameter of primary nozzle exit plane, in. 

DT 

d t 

diameter of throat, in. 

ETA 

T) 

efficiency, 100MN, percent 


g 

o 

acceleration due to gravity, 32.163 ft/sec 

GAMMA 

y 

specific weight of fluid, p(g/g c ), lbf/ft^ 



dimensional constant, 32.174 (ft-lbm)/(sec^)(lbf) 

KD 

K d 

friction loss coefficient for diffuser 

KP 

K P 

friction loss coefficient for primary nozzle 

KS 

K s 

friction loss coefficient for secondary inlet 

KT 

K t 

friction loss coefficient for throat 

M 

M 

flow ratio, Wg/W^ 

N 

N 

head ratio, (P D - P 2 )/( p 1 - P £>) 

NOMS 


number of values of M to be read as input 

NOP2 


number of values of Pg to be read as input 

NORS 


number of values of R to be read as input 

PI 

P 1 

primary total inlet pressure, psia 

P2 

P 2 

secondary total inlet pressure, psia 

P2REQD 

P 2REQD 

secondary inlet pressure required to avoid cavitation, psia 

PD 

P D 

outlet total pressure, psia 

PDREQD 

P DREQD 

total pressure to which jet pump is required to discharge, 
psia 

PV 

P v 

vapor pressure of secondary fluid, psia 
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FORTRAN Mathematical Definition 

variable symbol 


R 

STGMAL 


W1 

W2 

WT 

Subscripts: 

D 

d 

n 

P 

s 

T 

t 

1 

2 


R 

P 


V 

W 1 

W 2 

Wx 


area ratio, A n /A^ 
fluid density, lbm/ft^ 

jet pump cavitation prediction parameter at headrise drop- 
off, (P 2 - P v )/(rv2/2g) 
fluid velocity, ft/sec 
primary fluid weight flow, lbm/sec 
secondary fluid weight flow, lbm/sec 
total weight flow, + V^, lbm/sec 


discharge 

diffuser 

primary nozzle exit plane, jet pump 

primary nozzle 

secondary fluid inlet 

total 

throat 

primary fluid 
secondary fluid 


i 
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non 


APPENDIX B 


FORTRAN IV LISTINGS 


J FT PUMP L)ES IGN - I - 


D IM FN S If)\! R ( 30 ) * M < 30 ) , P2<20) 

REAL M.NTH.KP,KS,KT.KD,NN 

NAMFLIST/CAR01/KP,<S,KT,KD * GAMMA » S I GMAL t PV 
NAMEL I S T/ CAR 02/ P i,w2,NURS,RtNOMS*M*NQP2 ,P2 

HEAD l 5.CAR01 ) 1 

10 READ < 5.CARD2) 3 

W k I TE < 6, 100 I Pi, W2,PV, GAMMA • SI GMA L , KP • K$ ♦ KT . KD 4 

G AH 2G = GAMMA/64.326 
OU 1000 L=1.NDP2 

R I T E (6, 200 ) P2(L> T 


00 1000 J = 1 , N G M S 
C 

C CALCULATE FLOW RATES 

C 

M 1 = W2/M ( J ) 

WT = Mil + W2 
00 1000 l=l,NORS 
C 

C CALCULATE GEUMETRY, AN, DN, AT. ON. 

C 

AN2T = i.+KP - < i.+KS)*( (M(J)*R(I )/ll.-R(i ) ) )**2) 
IF (AN2T.LT.0. I GO TO 1000 


AN = < l*4.*Wl/64.348/GA02G)* SOK T ( AN2T / ( 144. *(P1-P2)/GA023) ) 2) 

ON = SJRTC AN/.7854) 21 

AT = AN/*< i I 
DT = SJRT( AT/. 78541 
C 

C CALCULATE heau ratio, efficiency, ano outlet pressure 

C 23 

M TH = NN( R ( I ) , M ( J ) ,KPtKS,KT,KO) 25 

IF (NTH.GE.O.) GO TO 999 

wRITF ( 6, 300) M(J),Wl«WT,R(I)« NTH 29 


GU TO 1000 
999 ETA = MU )* NTH 

PI) = ( mTH*P 1 + P 2 ( L ) ) / ( 1 • + NTH ) 

C 

C CHECK FOR CAVITATION 

C 

P2KEU0 = (SIGMAL* GAU2G /144.1* (( 1 44 . / ( 2 . *GA02G *32. 1 7 A ) 1**2)* 

K(((w? * Hill/ (AN *il. - RU)m**2) + PV 

IF <P?<L) .GT.P2REQDI GO TO 600 

WRITE (6, 400 )M< J), Wit WT,R< I ) . 3 N, 0 T , NTM , E T A, PD • P2 R EOO 41 

GO TO 1000 

600 d K I TF ( o, 500)M( J),W1, WT,R(I),DN,DT,NTH,ETA,PDf P2RE0D 45 

1000 CUNTINUE 
GO TO 10 
C 

C FORMAT STATEMENTS 

C 

100 FORMAT! 1H1.49X.21HJET PUMP DESIGN - I - // / /1 0 X , 2HP 1 ,7 X , 2 HH 2 , 7X • 

A? HP V, SX , 7 HGAMM A , 3 X . 7HS I GMA L , 5X , 29 KP , 8X , 2 HK S, 3 X , 2 HKT , 3 X , 2HKD' / 

RF13.l,F9.2, F9.1.2F1.0* l» 4F 10.3) 

200 FOKMaTI /// 9X.5HP2 = . F 5 . 2 , / /I OX . 1 HM , 8X , 2H Wl • 8 X . 2H WT , 8X , HR , 8 X . 

A2 HUN , 8X ,2 HOT «rt X , 3H N , 5X , 3HE TA • 7X , 2HP D , 5X , 6HP2RE 00 , 4X , 7HNC OR C) 

300 FORMAT ( F 1 3 . 1 . 2F 10 . 2 , F 9. 3 , 2 CX . F9. 3 ) 

400 FORMAT ( F13. 1, ?F 10 . 2 . F 9. 3 , F U. . 3, F 9. 3 , F 1 0. 3 , F9 . 3 , 2 F9 . 1 , 7 X , H C ) 

5GC F Ik MAT (F13. 1 . 2F 10 . 2 . F9. 3 , F 10. 3, = 9. 3 . F 1 0 . 3.F9.3 ,2F9.1 ,7X,2HNC) 

END 
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c 

C FUNCTION SlibPROGft A M FUR HEAD RATIO CALCULATION 

(. 

REAL FUNCTION NN(R*M,KP*KS,KT»ftD) 

REAL M. <P, K S* KT. KO 

TFRMi = Z.*R ♦ ((2.*RM2W*2)/(i,-(UN <1. * KT f KOI *R* *2 *( 1 . +M ) 
A** 2 

T FRM2 = 1(1* + KS)*R**2*M**2) /<1.-R)**2 
WN - /TE<Ml - TER M2 I /III. * KP) - TEkMl) 

< FT URN 
t NO 


C 

C it/ PUMP DESIGN - II - 

C 

DIMENSION ft 4 30 I t P2 ( 2G) • EM(30) 

R EAL M.NTH»KP»KS»KT » K3 , NN 

WAMFL IST/C ARDl/KP *<S*KT. K3 #GAMMA , SI SMALtP/ 

NAMFL IST/CAR02/P1, *2 • OT , NOP 2 * P2 ,NUkS. R 
READ I 5 , C ARD I ) 

10 READ (5.CARD2) 

wft ITE (6* LOU ) Pi ,W2 • OT • P Y » GAMMA * S I GMA L • • KS * KT , KO 
GAO 2G = i,AMMA/64. 326 
DO 1000 t =1,N3P2 
*kITF(6.200) P 2 ( L ) 

DO iOOU I - 1 * NO R S 

c 

C C Al CUL ATE GEUMETR i , AT* AN* UN* 

C 

AT= 0,7 d54*0T**2 

AN = AI *ft I I ) 

i)N = S(JRT( AN / O • 7B 54 ) 

0 

0 BEGIN APPROXIMATION PROCEDURE FOR PRIMARY F L 3 W RATE, 4 L • AND FLOW } AT I 3 , 

0 

ANUM = S-JkTl 144. MP1-P2IL) )/ GA02S ) 

FAC = IAM*64.348 * GA32S ) / 144* 

H NUM = ANiJM *E Al. 

DEN 1= SJ< T ( l .4-KP) 

FW1 = BNUM/UtNl 
EMI 1 ) = w2/E w 1 

(. 

L ENTER LOOP FOR CALULAT10N UE FLUw RATIO. 

c 

DU 20 Ou J = 1 * 1 5 

8 D E N = I l .*KP >-< 1 . + KS 

IF (HOEN.LT .0. > GO TO dOO 

CUfcN=SOKl I riUt-N I 

«Wl= riNtJM /COhN 

E M ( 0 ♦ i I =W 2 / w rt 1 

DEI Ms AttSl + J) ) 

ZhlHljkV = Del 4/EMI JFl) 

IF ( PFR.JhN .LE .0.0005) M=E.MI J*l ) 

(FI PEOEV .Lb .u.0005) GU TQ 2 > 3 
2000 CONTINUE 
Gil Tl) 7 50 
C 

c cai.uj.atc Primary elo* rau, mead ratio, efficiency, and ojilet pussj^e. 

r. 

2bu El)Fn= l. + RP - ( 1 • FKS ) * I M# K I 1 ) / ( I.-RI I) ) ) ** 2 
I E I EUEN.L 1.0.) i 3 TO 850 
E uEN = SJR r ( EJrN ) 

Wla 8 4UM/F0EN 



r> r> r> 


MTN = mM k ( I )*M.K 3 »K5.KT,K3) 45 

iH mth.o* .ii . ) so r;j 999 

«RITP (6.300k \i H .DN.Wl .NTH 49 

GU Tfl X 0(1 0 
999 ETA a M*NTH 

P 0= < NT H^P L + t>2(U)/( i.fNH) 

C 

C CHECK FiK LAVITATIJN 

c 

P?REuh> = (SIGNAL* GAQ2G /144.)* (< l 44 . / ( 2 . *0 A02 G *3 2 . I 7 V N ** 2 > * 

X(!(w2 * X ( I * )/ (AM *<i. - RU*m**2l * PV 
IF (P2RFJD.LT.P2) GU TO 45C 

M R I T t !t?,4lOJ 4 ( I | .ON.wl .N1H. M,=TA,P3,P2REwD 53 

(>) T) ID 3 0 

450 * k l T E (6,400) R ( t ) . DM . w I . N TH. M . E T 4 , P 3 . P2 RE UO 61 

Si) TH 103 0 

750 wKiTtf 6.751) R( I ) 64 

i.n tr 109 0 

600 rik HE (bfdOll k(l) 67 

G9 n 1090 

860 *R I T E <6,651 ) R (I ) 73 

13 Ol CONTI.Mjfc 
GO Tn 10 

c 

C FORMAT ST AT fc ME MTS 

L 

100 FORMAT! iHl .49X. 2HJ8T PUMP 3=5ISN- II -/ / I / 1 3 X , 2 H P 1 , 7 X , 2 -U l , 7 < , 

A2HGT. 7X .2HPV.5X ,7iGAMMA , j X , 7 FtS I 3 MA l .5Kf2H<P,&X.2HKS.3(.2H<T.3X, 

62FK0//F 13.1. F9. 2.-9. 3. F 9. 1 ,2F LO. 1 . 4F 13. 3) 

200 FORMAT! // /9X.5HP2 = , F 5 . 2/ / 1 2 X ,1 MR . 6X , 2 HJ N* 3 X , 2 H Wl • 9 X , 3 H M ,6X. 

AlHM.6X,3MFTA.6X.2MP3,5X,6HP2kEOJ.4X.7H.\lC JR C / / » 

3 00 FORMAT ( F 1 5 . 3. F 8 . 3. f 9* 2 . F U • 3 ) 

A 00 FORMAT (FlD.3*r8.3,-9.2» r ll. 5, P 6, ? *F3*i*F7.l»F9.1 , 9 < * 2 H M3 ) 

410 FUMAT ! F 15. 3,= 6. 3,- 9. 2.-LI . J . = 8. 2 .F) . 3 ,-7. I ,F9. 1 .9< . IhC) 

751 FORM A T(Fl5.3,10X, 30HPR 3GRA M UNABLE TO CVLCJLATE Wl) 

R31 FORMAT! hi 5.3. 10X. 15HM6G SORT IM HH 1 * 

651 FORMAT! F15.3.10X. 14HNE^ SORT IN Wl > 

END 


r. 

c FJNCTlON SUBPROGRAM FOR FiF AO RATIO 34L3 JL AT I ON 

C 

RFAL FUNCTION NN( R * M,KP,KS*KT, KD) 

«FAl M, <P, KS. KT, X9 

IFRM1 = 2 • *k + ! (2.*R**2*M**2) /< i.-R) )-( L.HCT fKO) *R** 2*( i .PM) **2 
term? = ( ! l.*KS)*R**2*M**2)/U.-R)**2 
N N = iTFKMl - TERM2) /! < 1 . P KP ) - TfcKMl ) 

RFTUR M 
END 


JET FLMP C85IGN - III - 

C I N EK 5 1 ON M30), R(3C),P2(30) 

REAL f,MFfKP*KSfKT f KDfNN 

NAFEL I ST/CARD l/KP f KS,KT, KD, GAMMA , SIGMAL f PV 
(NAf'EL IST/CARC2/P 1»PCREQD,W2.N0MS,M t NORS, R 

REAC ( 5 » C £ R C 1 ) 1 

C-AC2G = CAPPA/64.326 
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ooo non non nnon on o 


1C RE AC ( 5 * CA RC 2 ) 3 

WR I T E ( 6 » ICC ) P1,PCRECD,W2,KP,KS,KT,K0»GAMMA .SIGMAL.PV A 

CC 1CC0 1 * 1 »NOMS 

WRITE(6,15C) MO) 7 

Wl = Vi2/ MI) 

CC 1CC0 J* 1 , NORS 

TEM = ( (M t)*R(J )/( l.-RI J) ) )**2) * (l.+KS) 

CALCULATE HE AO RATIO AND EFFICIENCY 

NTF* NNI R ( J ) , Mil), KP.KS*K T »KD ) 16 

IFINTH.GE.C.) GO TO 200 
GC TC 10C0 

2CC ETA = NTH * Mil) 

BEGIN ITERATION LOOP ON P2. FIRST ESTIMATE OF P2 CALCULATED FROM 
CEFIMTICN CF HE AC RATIO. 

P2 ( 1 )* PCRECC - NTH ♦ (P1-P0RE0D) 

DC 25C 1*1,15 

C * (144. AM) / ( 32. 174*GA02C * 2.) 

Cl * 1. i KP - TEM 

IF (C1.LI.C.) GO TO 750 

C2* (PI - P 2(D) * 144./GA02G 

»1* C*SCRT(Ci/C2) 30 

P2RECC * (SIGMAL* GA02G /144.J* (( 144. /(2.* GA02G*32 . 174 ) ) **2 ) * 

X ( I (W2 * F ( J ) ) / ( AN1*( 1. - R ( J > ) J )**2 ) ♦ PV 

CELP* P2RECC - P2(L) 

PEPCEV * A8S< 100.*DELP/P2REQC ) 

IF ( FERDEV .GT. .05) GO TO 24C 

WHEN FERCENT OEVIATICN IS LESS THAN OR EQUAL TC .C5, END THE LOOP 

P2F * P2FECC 
AN = AN 1 
GC TC 3CC 

24C P2(L«1I* P2REOO 
25C CONTINUE 
GC TC 8CC 

CALCLLATE CUTLET PRESSURE ANC GEOMETRY, AT, DT, ANC DN. 


3CC PC* (NTH * PI ♦ P2F )/( l.+NTH) 

AT * AN/FIJ) 

CT = SORT (AT/ .7854 ) 45 

CN* SCRT ( AN/ . 7854 ) 46 

WRITE (6,110) RtJ), NTH, ETA, P2F , PO, DT, DN, W1 47 

IFIPC.LT. PCRECD) GO TO 700 
GC TC 1CC0 

7CC WRITE(6,12C) 53 

GC TC 1CCC 

75C WRITE (6,135) 55 

GC TC ICCO 

8CC WR1TE(6»13C I PERCEV 57 

10CC CONTINUE 
GC TC 1C 


FCRMAT STATEMENTS 

ICC FCRMAT ( 1 H 1 47X, 22HJET PUMP CESIGN- III -////13X ,2HP1 ,5X ,6HPDRE0D, 
A5X,2FW2,7X,2HKP»i8X,2FKS*8X,2HKT,8X,2HKD,eX,7HGAMMA ,3X,7HSIGMAL , 
e5X, 2HPV// F 16. 1 , 2F9.1, 4F 1C. 3 , 2F 10. 1 , F9.1 ) 

11C FCPMAT (F17.3,F10.3,F9.3,F10.1,F11.1,2F11.3,F10.2) 

12C FORMAT (1H+ 92X, 22HPD IS LESS THAN PORECD ) 

13C FCRMAT <1M 24X^ 1 2HPERDEV (15) *, F 1C . 8 , t X , 3 1HC0NVERGENCE ON P2 D1 
AC NCT OCCUR ) 

135 FCRMAT ( 1H+ 30X, 31HC0NVERGENCE ON P2 DIC NCT OCCUR I 
15C FCRMAT (/// 12X,4HM = , F6. 3/ / 14X , 1HR , 9X , 1HN,7X, 3HE TA , 5X , 7HP2 REU 
AC, 6X, 2HFC, 8X, 2HCT , 9X, 2H0N, 9X.2HWI //) 

ENC 
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c 

C FACTION SI E PROG RAM FOR HEAD RATIO CALCULATION 

C 

REAL FUNCTION NN < R , M t KP , KS # KT, KD ) 

REAL K, K P , K S, KT » KD 

TERM = 2.*R + { (2.*R**2*M**2)/( l.-R ) )-( 1 . +K T+KD > *R **2* ( l.+M)**2 
TEPR2 = l(l.+KS)*R**2*M**2>/U.-R)**2 
NN = (TERM - TERM2) /((l. + KP ) - TERM 1 ) 

RETURN 

ENC 


C 

C J FT PUMP OES IGN - IV - 

C 

DIMENSION R( 301 * MOO) 

REAL M.NTH.KP, KS.KT.KD.NN 

N AM EL IST/CARDl/KP.KS.KT* KD » G A M M A * SI GMAL , PV 
NAM EL I ST/ GAR D 2 /N 0 R S , R » NO MS ,M,P 2 ,PD,W 2 


READ (5.CAR01) 1 

10 READ < 5* C AR 02 ) 3 

WRITE (6,100) P2, ? D,W2tPV * GAMMA * S I 3 MAL , KP , KS , KT , K 0 A 

GAU2G = GAMMA/64.326 
00 1000 J = 1. NOR S 

W R ITE (6, 200 I R< J) 7 

no 1000 I 3 1 , NO M s 

Ml* W?/M< I ) 

NTH* NN(R( J).M( I > . KP . K S , K T , KD ) 13 

IF (NTH.GE.O.) GO TO 999 

WRITE (6,300) M(I).W1,NTH 17 


GO TO 1000 
999 ETA* M ( I1*NTH 

Pi = (POP2) / NTH ♦ PD 
C 

c calculate geometry, an. at, on, ot. 

c 

AN2T * 1. ♦ KP - ( i.+KS)*( ( M(I)*RI J) / ( 1 . -R ( J )) ) * * 2 ) 

IE ( 4N2T.LT.0. ) SO TO 1000 

AN=< l44.*wl/64.348/GA02G)*S0RT (AN2T/(144. *( P1-P2 ) / GA02 G > ) 28 

AT* AN / R( J ) 

OT* SORT( AT/. 7854) 3D 

0N= SO 3 T ( AN/.7854I 
C 

C CHECK EOT CAVITATION. 

C 31 

P2REGD * ( S I GM A L ♦ G A202G/ 1 44. ) * (( 144. /(2.* GA02S*32. 1 74 ) ) ** 2 ) * 

X(((W2 * R ( J ) )/ (AN *(i. - K(J))))**2> + PV 

I F( P2 .GT. P2RE00) GO TO 550 


WRITE 

< 6 . 399 ) 

M ( I 1 » 

Wl, 

NTH. 

ETA, 

PI . 

AN. 

AT, 

DT, 

ON 

3 S 

GO TO 
550 WRITE 

loan 

( b, 400 ) 

M ( 1 ) , 

wit 

NTH. 

ETA, 

PI , 

AN. 

AT, 

DT, 

DN 

39 


1000 CONTINUE 
GO TO 10 


f. 

C FORMAT STATEMENTS 


100 FORMAT < 1HI 49X.21HJET PUMP DESIGN- IV -// / / 10 X, 2 HP2 , 7 X , 2 HPD , 
A7X.2HW2.7X.2HPV. 5X, 7HGAMMA .3X.7HSIGMAL , 5X . 2 HKP , 3 X , 2 HKS . 3X, 

82 HK T . RX . 2 HKD //F13.2.3F9.2.2F10.2.4F10.3) 

200 FORMAT (///9X.4HR = , F 5. 3 . //I OX , 1 M M. 8X , 2H Wi , 7 X , 3H N ♦7X.3HETA.7X, 
A2HP1.8X.2HAN.10X.THNC UR C.5X.2HAT.9X.2HDT.9 X .2HDN// ) 

300 FORMAT ( F 1 3 . 2, F 9 . 2 , F 10. 3 ) 

399 FORMAT ( E l 3 . 2, F 9 . 2, 2F 1C . 3 , F 10 . I , F l C. 3 . 1 0 X , i HC , 7 X. FS . 3 . 2 F i l . 3 I 


400 FORMAT* ( F 1 3 . 2. F 9. 2. 2F 1 C. 3 . FI 0. 1 , F 1 0 . 3 . 10 X . 2H NC , 6X ,F 6, 3 .2 F 1 1 . 3 ) 

END 
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c 

C FUNCTION SUBPROGRAM FUK HEAD RATIO CALCULATION 

C 

REAL FUNCTION MN ( I . M . KP • KS * KT • KJ ) 

REAL M « <P* KS* KT* KO 

TERM1 = 2 .* k «- ( (2.*R**2*M**2) / ( 1 .“R ) I - ( 1 • + KT + KD ) *R**2 *(L.*M|**2 

TERM? = < ( 1.fKS)*R**2*M**2)/<1.-R)**2 

si N = ( TER M 1 - TERM2) /HI. * KP) - TERM1 ) 

i FTURN 

END 


C 

C JET PUMP OESIGN - V - 

C 

DIMENSION M ( 3C ) 

< FAL M * 4 *KP *KS * KT « KO • NN 

NAMFL IST/CAKDI/ Kf> .KS .KT . KD *G A MMA .SIGMAL.PV 
NAMfcL l ST/CAkl)?/ P l . P 2 * U N . 0 T . R *NOMS *M 
R F AO ( 5. C ARD 1 ) 

10 READ (5.CARD2I 

• RITE I 6* 5QC t Pl*P2.k*ON*OT*KP.KS*KT*KO, GAMMA .SIGMAL. P/ 

GA02G = UAMMA/64.326 
AN = .7854 * ON **2 
DO 1000 1=1* NO M S 

N = NN (R,M(i), KP*KS*KT.KU) 

F T A = N * M ( I ) 

PD = ( N*P 1 ♦ P? ) / ( 1 • + \l ) 

C 

C CALCULATE PRIMARY FLOW R \ TE • Wl 

C 

TEM = 1. + KP - II. * KS) * (R*MU)/U. - R>)**2 
IF(TEM .LF:. 0.) GO TO 600 

Wl = I GA02G* 2.* 32. 174*AN/144. ) * SUR T( ( I P i-P 2 ) * 1 44 . i GA02 3 ) t T E M ) 

H? = wl * M ( l ) 

C 

L CHECK FUR CAVITATIUN 

C 

P2REQ0 = (SIGMAL* GAU2G /144.)* (I 1 44. / ( 2. *GA02G *32. 1 7 4 I ) **2 ) * 
X(((W? * R )/ IAN * ( 1. - k )l 1**2) F PV 
IF (P2REJD.LT. P2) GO TO 700 
•RITE (6.520) M( l ) . N.fcTA *PQ *fcl ,M2 , P2RE03 
GO TO 1000 


600 •RITE (6. 530) 

GO TO 1030 

700 WRITE (6.510) M (I ) . N . E TA , PD . M 1 * *2 • P2 R EGD 
1000 CONTINUE 
GO TO 10 


C 

C FORMAT STATEMENTS 


600 FORMAT (1H1 47X. 21HJET PUMP DESIGN - V - / /V / 4 X . 2HP 1 , 5 X , 2 H P2 . 6X * 
AlHK. 6X.2HDN.5X. 2r|3T, 6K . 2 HK P .6X . 2 i* S. 6 X . 2 HKT * 6 X , 2 HKD . 4X . 7 HG A MMA 
B . 3X * 7HSIGMAL »5X. 2 HP V / / F8 . 1 . F 6. 1 *F 7 . 3 . F8 . 3 , F7 . 3 . 4F 8 . 3 . F 8 • 1 , 

C2F 10. I // /13X. IHM. 9X. IHN . 8 X . 3HE T A . 8X , 2 H P3 . 9X . 2 H Wl * 9 X. 2 * 9 X * 6 HP2 

ORFOO // ) 

510 FORMAT ( 10X . F 6* 2 . 2F 1 0. 3 . E1G.2. 2F11.2 .F12.1) 

520 FORMAT ( 10 X . F 6. 2 . 2F 1 0 . 3 . F10.2. 2F11.2 »F12.1,2H C) 


530 FORMAT ( 1 OX • 20HNE G A T I VE SUUARE ROOT 


) 


FNO 


9 


15 


19 

22 

2V 
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RJiViCTUJN SUbFRUURAM F UR HEAD RAT 13 C A LuULAT I UN 


REAl FUNCTION MN<R*M*KP*K$ # KT «KL)I 
REAL M. <P* KS* KT * KD 

TERMl = 2.*R + ((2.*R**2*M**2)/(l.-Rn-(l.*KTM«)l*R**2*<l.«-M»**2 
TERM2 = ( ( l.+KS)*R**2*M**2l/(l.-R)**2 
N = ( TER M 1 - TERM2) /(Cl. ♦ KPJ - TERMl ) 

RETURN 

END 
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